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1.  INTRODUCTION 

Most  modern  structural  adhesives  are  designed  for  high  tough- 
ness and  therefore  exhibit  highly  non-linear  inelastic  character- 
istics. Hence,  a non-linear  program  is  essential  for  predicting 
structural  performance  of  a bonded  system. 

Due  to  the  complex  visco-elastic-plastic  nature  of  the  IAL 
even  at  low  stress  levels,  rigorous  analytical  formulation  of 
its  mechanical  behavior  up  to  failure  is  impossible,  especially 
when  time-  and  environmental  factors  are  taken  into  account. 

Thus,  comprehensive  empirical  information  is  necessary  on  the 
inelastic  stress-strain  relationship  of  the  adhesive  as  a func- 
tion of  time,  temperature  and  such  factors  as  moisture  content. 

Bulk  data  on  the  adhesive  may  suffice  so  long  as  stress 
distribution  and  cohesive  failure  within  the  bond  are  concerned. 
However,  when  failure  initiates  or  propagates  close  to  the 
adhesive-adherend  interfacial  zone,  in-situ  characteristics 
must  be  taken  into  account  and  correlated  with  the  bulk  properties, 
Accordingly,  a series  of  tests  on  the  in-situ  mechanical  behavior 
of  an  adhesive,  for  different  loading  modes,  geometries  and 
strain  rates,  was  conducted  with  a view  to  correlating  the  find- 
ings with  the  behavior  of  the  bulk  adhesive  under  uniaxial  load. 
This  goal  was  partly  achieved  and  results  are  described  in  the 
present  report  which  deals  with  three  main  subjects: 

(a)  The  mechanical  characteristics  of  the  bulk  epoxy 
used  in  the  adhesive  layer. 

(b)  The  in-situ  mechanical  characteristics  of  the 
adhesive  layer. 

(c)  The  two-dimensional  stress  distribution  (analytically 
solved)  within  the  IAL  at  the  non-linear  range  of 
the  adhesive. 

Preliminary  tests  on  the  in-situ  hygrothermal  behavior  of  the 
adhesive  will  be  reported  briefly  in  conjunction  with  a future 
research  program. 

The  following  is  a list  of  reports  related  directly  or 
indirectly  to  the  output  of  the  three-year  research  activity 
sponsored  by  the  U.S.  Army;  they  have  either  been  published 
or  submitted  for  publication. 
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2 . ADHESIVE  BULK  CHARACTERI ST.,  JS 

2.1  Introduction 


Bulk  epoxy  characteristics  under  different  loading  modes  and 
testing  conditions  are  currently  obtained  using  standard  experi- 
mental facilities  such  as  the  Instron  tester.  A ductile  type 
of  epoxy-versamid  system,  which  was  previously  investigated  [1-4], 
was  selected  to  serve  as  representative  for  structurally  tough 
adhesives  throughout  the  present  investigation.  Most  of  the  data 
is  available  for  uniaxial  tension.  The  purpose  of  the  present 
test  series  was  to  relate  this  data  to  other  loading  modes  such 
as  shear.  This  relation  will  enable  the  prediction  of  the 
general  stress-strain  behavior  of  bulk  epoxy  from  the  elastic 
stage  through  the  inelastic  one  up  to  failure,  under  a combined 
state  of  stress. 

2.2  Specimen  Fabrication 

The  material  tested  consisted  of  Shell  epoxy  epon  815  and 
General  Mills  versamid  V-140  in  the  ratio  of  70:30.  The  mixture 
was  cured  for  24  hours  at  room  temperature  followed  by  post- 
curing of  6 hours  at  80°C. 

For  the  tensile  test,  dog-bone  specimens  were  cut  from  the 
cast  plates  according  to  ASTM  D638-64T.  For  the  shear  test, 
thin-walled  tubes  were  prepared  and  machined  to  have  thinner 
thickness  along  the  gauge  length.  Special  metal  grips  were  pre- 
pared and  bonded  to  the  specimens  as  shown  in  Fig.  4. 

2.3  Uniaxial  Tension  Loading 

The  specimens  were  loaded  m the  Instron  tester  under  uni- 
axial_tension  at  various  strain  rates  ranging  from  l5xl0_,*min-x 
38x10  z at  room  temperature  (about  22 “C)  A typical  stress- 
strain  curve  is  shown  in  Figure  1.  The  general  trend  for  all 
curves  shows  a short  linear  portion  beyond  which  a non-linear 
curve  terminates  in  a stress  plateau. 

Other  specimens  were  similarly  loaded  in  uniaxial  tension, 
and  transverse  strain  was  recorded  simultaneously  with  longi- 
tudinal strain  by  means  of  a special  extensometer . The  ratio  of 
transverse  to  longitudinal  strains  provides  the  Poisson  ratio, 
v,  which  is  plotted  as  a function  of  uniaxial  strain  in  Figure  1. 
The  strain-rate  effect  can  be  evaluated  from  the  family  of  stress- 
strain  curves  obtained  under  different  strain  rates  (Fig.  2). 

The  major  effect  is  on  the  maximum  stress  plateau,  S^,  which  can 
be  defined  as  yield  stress,  and  shows  an  increase  miineanty 
with  log  strain-rates  (1]  (Fiq.  3).  The  effect  on  initial 
Young's  modulus  and  proportional  limit  is  less  pronounced  (Fig. 3). 


Tubular  specimens  were  mounted  on  an  Instron  tester  and  loaded 
in  torsion.  A device  was  designed  to  measure  change  of  torsional 
angle,  using  the  regular  Instron  extensometer  (Fig.  4).  The 
overall  torsional  deformation  of  the  specimens  consists  of  an 
unknown  displacement  at  the  grip  region  which  is  difficult  to 
estimate.  In  order  to  obtain  the  net  shear  deformation  along  the 
uniform  gauge  length,  the  following  procedure  was  carried  out: 
Tubular  specimens  of  different  lengths,  ranging  from  40  to  110  mm, 
were  loaded  in  torsion  up  to  failure.  Torsional  displacement  at 
each  torsion  moment  level  was  plotted  as  a function  of  specimen 
length.  These  plots  show  a series  of  lines  (Fig.  5)  which  can 
be  extrapolated  to  intersect  the  zero  length  axis  at  different 
torsional  displacement  levels.  The  extrapolated  shear  strain 
at  zero  length,  as  a function  of  the  respective  shear  stress 
(calculated  from  the  respective  torsional  moment)  provides  the 
correction  curve  (Fig.  6) . By  subtraction  of  the  correction 
curve  from  the  overall  shear  stress-strain  curve,  the  "true"  shear 
stress-strain  relationship  can  be  obtained  as  shown  in  Figure  6. 

2.5  Effective  Stress-Strain  Relationship 

An  effective  stress-strain  relationship  (Fig.  7,  solid  curve) 
was  derived  from  the  uniaxial-tensile  stress  (Fig.  1)  for  the 
representative  epoxy  resin.  The  effective  stresses  and  strains 
are  related  to  their  respective  stress  components  by  using  the 
von  Mises1  deviatoric-er.ergy  yield  criterion  [5]  as  follows: 


s = C [a  -a  )2  + (o  -o  )2  + (o  -a  ) 2 
x y y z z x 
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In  the  case  of  uniaxial  tension  test  data,  the  effective 
stress-strain  curve  is  identical  to  the  respective  uniaxial 
tensile  stress-strain  curve  of  Figure  1.  In  order  to  calculate 
the  effective  stress-strain  relation  from  the  shear  tests  data 
(Fig.  6)  the  following  relations  are  used: 


= /3 


xz 


e = 


/3  r xz 


2 1+v 


(3) 


It  has  to  be  noted  that  in  order  to  compare  shear  data  with 
respective  tensile  data,  one  must  use  the  same  effective  strain 
rate , i . e . , 


"XZ  a+0)  /3  "x 


xz 


/3 


(4) 


• • 

where  yxz  and  sx  are  the  shear  and  tensile  strain  rates  respec- 
tively. 

The  plots  of  shear  and  tensile  tests  data  in  terms  of  the 
effective  stress-strain  relationship  in  Figure  7 exhibit  good 
agreement.  This  supports  the  effective  stress-strain  concept 
used  in  the  present  study. 

The  analytical  solution  for  inelastic  stress  distribution 
in  the  adhesive  within  the  doubler  model  necessitates  the  deriva- 
tion of  the  adhesive  moduli  as  well  as  Poisson  Ratio  at  the  dif- 
ferent strain  levels. 

The  variations  of  effective  tangent  (E^)  and  secant  (Es) 
moduli,  as  well  as  effective  Poisson  ratio,  v,  as  functions  of 
effective  strain,  are  shown  in  Figure  8.  Three  regions  may  be 
distinguished  along  the  stress-strain  curve,  namely: 

1.  an  almost  linear  region  up  to  the  proportional  limit 
(Sy)  through  which  v is  almost  constant; 

2.  a non-linear,  probably  viscoelastic  range,  characterized 
by  an  increase  in  v,  up  to  the  "yield  plateau"  level 

(Sf ) ; 

3.  a macro-plastic  flow  range  beyond  which  stresses  are 
almost  constant  and  even  tend  to  drop  at  higher  strain 
levels.  v is  almost  invariant  at  this  range.  By  sub- 
straction  of  the  elastic  components  of  axial  and  trans- 
verse strains,  the  "non-elastic  Poisson  ratio"  may  be 
derived.  At  the  "yield  plateau"  its  value  tends  to 
approach  0.5,  which  is  the  theoretical  limit  based  on 
the  assumption  of  incompressible  flow. 


— 
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All  the  characteristics  described  in  Figures  7 and  8 are 
functions  of  loading  strain  rate  and  temperature.  The  present 
study,  however,  will  be  confined  to  the  reference  strain  rate 
and  temperature  specified  in  Figures  1 and  7 which  are  considered 
to  be  typical  and  representative  for  similar  "non-linear"  ad- 
hesive materials. 

2.6  Bulk  Characteristics  Compared  with  Adhesive  In-Situ  Data 

In  a previous  report  and  subsequent  publication  [6]  an  inves- 
tigation of  shear  stress-strain  relationship  of  a similar  epoxy 
adhesive  in-situ  was  reported.  The  study  was  concerned  with  the 
effect  of  adhesive  layer  thickness  on  its  mechanical  characteris- 
tics. The  main  objective  was  to  compare  the  behavior  of  the 
adhesive  layer  in-situ  with  its  bulk  reference. 

Stress-Strain  Behavior 

Shear  stress-strain  curves  for  adhesives  of  different  thick- 
ness compared  with  the  bulk  references  are  shown  in  Figures 
6 and  9,  The  similarity  of  bulk  to  in-situ  behavior  is  evi- 
dent at  the  linear  elastic  ranqe  up  to  the  initiation  of  the 
yield  plateau  level. 

Adhesive  Shear  Elastic  Modulus 

The  effect  of  in-situ  adhesive  layer  thickness  on  its  shear 
modulus  is  shown  in  Figure  10.  The  level  of  bulk  modulus  seems 
to  be  approximately  the  same  as  the  thick  adhesive  layer.  It 
may  be  concluded  that  no  significant  variation  in  modulus  is 
apparent  for  adhesive  layers  of  thicknesses  above  0.3mm. 

Adhesive  Shear  Strength 

The  effect  of  in-situ  adhesive  layer  thickness  on  its  shear 
strength  is  shown  in  Figure  11.  The  shear  strength  of  the  bulk 
epoxy  seems  to  be  on  a slightly  lower  level  compared  with  that 
of  the  adhesive  in-situ.  It  may  be  concluded  that  for  epoxy 
adhesive  thickness  above  0.15mm,  shear  strength  is  almost  in- 
variant with  adhesive  thickness.  Hence,  bulk  shear  data  can 
be  used  for  a first  estimation  of  bond  strength  as  long  as 
failure  is  of  the  cohesive  type  within  the  adhesive  layer. 
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3.  MECHANICAL  CHARACTERIZATION  OF  ADHESIVE  LAYER  IN-SITU 
UNDER  COMBINED  LOAD 

3.1  Introduction 


Knowledge  of  the  mechanical  properties  of  adhesives  in  a 
combined  state  of  stress  is  essential  in  engineering  applica- 
tions of  structural  bonded  joints  in  which  the  adhesive  layer 
undergoes  a complex  state  of  stress  even  under  simple  loading. 

Unfortunately,  experimental  data  in  this  content  is  very 
scanty.  Early  studies  of  the  single-lap  joint  model  [7]  tended 
to  oversimplify  the  problem  by  assuming  pure  shear;  subsequently, 
the  more  sophisticated  closed-form  solution  of  Goland  and  Reissner 
[8]  envisaged  a lateral  normal  stress  component,  the  so-called 
"peel  stress",  while  more  recent  experimental  work  deals  mostly 
with  cases  of  pure  shear  [6]  or  tension  [9] 

The  investigation  reported  here  was  conducted  using  a torsion- 
tension  apparatus  to  measure  the  mechanical  properties  of  thin 
adhesive  layers  in-situ.  Some  preliminary  results  are  presented, 
indicating  the  suitability  of  the  method  for  the  purpose  in 
question . 

3.2  Testing  Apparatus  (Figs.  12a, b) 

The  test  specimen  consists  of  a pair  of  2025  aluminum 
flanged  cylinders  (adherends)  a.  The  loading  device  consists 
of  a torsion  arm,  b,  operated  from  any  conventional  loading 
tester  (e.g.  Instron) , and  a spring  operated  axial  loading 
post,  c,  capable  of  applying  tension  or  compression 

3.3  Procedure 


The  test  adhesive  was  composed  of  Shell  epoxy  resin  815  and 
General  Mills  versamid  V-140  in  the  ratio  of  70:30.  The  adherends 
were  aligned  on  a special  fixture  with  a controlled  gap  of  1mm, 
into  which  the  resin  was  poured.  The  joint  was  cured  for  24 
hours  at  room  temperature,  plus  6 additional  hours  at  60°C.  The 
specimens  were  then  tested  on  a conventional  Instron  machine  at 
different  constant  shear-strain  rates  and  under  different  axial 
load  levels. 

For  calibration,  shear  displacements  of  specimens  with  "zero 
adhesive  thickness"  were  recorded  and  deducted  from  the  overall 
shear  displacement  of  the  bonded  specimens.  The  shear  strains 
and  stresses  are  computed  from  recorded  moment-displacement  curves. 

3 . 4  Results 


Typical  shear  stress-strain  curves  in  Figure  13  indicate 
that  normal  uniaxial  load  has  a clear-cut  effect  on  the  inelastic 
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behavior  of  the  adhesive.  Axial  tension  tends  to  reduce  the 
initial  modulus,  shear  strength  and  ultimate  strain,  while  axial 
compression  tends  to  increase  them.  However,  because  of  the 
lower  in-situ  tensile  strength  of  the  adhesive-adherend  inter- 
face and  its  brittle  mode  of  failure  under  tension  load,  the 
range  of  applicability  is  limited.  Stiffness  and  level  of 
yield  plateau  seem  to  be  lower  with  the  reduction  in  strain  rate. 
In-situ  stress-strain  relationship  shows  fair  agreement  with 
reference  bulk  data.  Failure  under  combined  load  may  be  des- 
cribed by  means  of  the  apparent  ultimate-stress  combination 
plotted  in  the  shear-normal  axial  stress  plane  (Fig.  14) . The 
failure  envelope  thus  obtained  also  shows  reduction  of  the  shear 
strength  with  increasing  tensile  load,  and  increase  under  com- 
pression . 

It  may  be  concluded  that  the  above  methodology  for  adhesive 
in-situ  characterization  under  combined  loading  can  provide  use- 
ful data  for  failure  analysis  of  bonded  systems. 
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4.  THE  EFFECT  OF  HYGROTHERMAL  CONDITIONS  ON  MECHANICAL 
BEHAVIOR  OF  THE  ADHESIVE  LAYER. 


The  above  topic  was  scheduled  to  be  initiated  only  at  an 
advanced  stage  of  the  research  due  to  its  complexity  and  the  need 
for  adequate  experimental  and  analytical  basic  substrate.  It  is 
aimed  to  achieve  a more  realistic  picture  of  IAL  behavior  under 
external  conditions,  and  also  to  provide  some  information  on  the 
long-term  durability  of  structural  adhesive  joints. 

At  the  present  preliminary  stage,  the  effect  of  hot  water 
immersion  on  the  adhesive  layer  was  investigated  and  compared 
with  an  oven  dried  counterpart.  Preliminary  results  (Fiq. 

15)  indicate  a significant  change  in  the  shear  stress-strain 
characteristics.  The  presoaked  adhesive  specimen  showed  lower 
stiffness,  strength  and  ultimate  strain  compared  to  its  pre- 
viously dried  reference. 

The  next  step  will  be  a more  comprehensive  investigation  of 
adhesive  specimens  exposed  to  different  hygrothermal  conditions 
and  tested  in  shear  and  tension  in  wet  and  dry  conditions. 

The  long  range  scheme  is  to  measure  hygrothermal  effects  on 
interlaminar  stress  and  strain  distributions  within  a doubler 
specimen  by  utilizing  the  method  described  in  the  Final  Techni- 
cal Report  of  the  second  year.  The  experimental  results  will  be 
compared  with  the  analytical  FEM  solution  and  will  take  into 
account  the  hygrothermal  effects  on  the  adhesive  stress-strain 
relationship . 
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5 . INTERLAMINAR  STRESS  DISTRIBUTION  WITHIN  ADHESIVE 

LAYER  AT  THE  NON-LINEAR  RANGE. 

5 . 1 Introduction 

Numerous  publications  are  available  on  the  analytical  solu- 
tion for  stress  distribution  within  the  adhesive  layer  of  a 
structural  bonded  joint  model.  Most  of  the  earlier  works  [8,10- 
13]  are  based  on  certain  assumptions  such  as  uniform  shear  and 
normal  stress  distribution  through  the  adhesive  thickness  which 
enable  the  derivation  of  a closed  form  solution.  These  simpli- 
fications, which  provide  only  the  average  stress  data,  lead  also 
to  certain  inconsistencies  in  the  aquilibrium  equations  and 
violate  boundary  conditions  at  the  edges. 

More  recent  works  deal  with  the  analytical  complexity  of 
the  problem  by  applying  numerical  methods  such  as  finite  dif- 
ference [14]  and  finite  element  [15-18]  and  provide  the  two- 
dimensional  stress  distribution  within  the  adhesive  layer  for 
various  bonded  joint  models  of  isotropic  [15,16]  and  orthotropic 
[17,18]  adherends.  These  solutions  are  in  agreement  with  equili- 
brium and  boundary  conditions  of  the  structural  model  but  are 
limited  to  the  elastic-linear  stress-strain  relationship. 

Unfortunately,  the  polymeric  adhesive  layer  commonly  used 
is  characterized  by  a non-linear  stress-strain  behavior,  even 
at  a relatively  low  stress  level.  This  non-linearity  is  more 
pronounced  at  high  stresses,  reaching  a stress  plateau  which 
may  be  defined  as  "macro-plastic"  yielding  (see  Fig.  7).  Another 
approach  to  represent  this  mode  of  behavior  is  by  a simpli- 
fied elastic-plastic  model  [ 1 8 - 2 u j wmch  may  provide  an  approxi- 
mate upper  bound  for  the  real  solution.  Such  a solution  is 
insufficient  considering  the  fact  that  under  service  load,  the 
material  is  not  allowed  to  approach  this  "yield  point"  and  that 
at  the  critical  region  of  the  adhesive,  the  stress  level  is 
mainly  at  the  non-linear  range,  beyond  the  elastic  limit  but 
below  the  yield  plateau. 

The  viscoelastic-plastic  nature  of  the  polymeric  adhesive, 
which  is  reflected  in  its  sensitivity  to  temperature  and  its 
dependence  on  time  [3,4],  also  contributes  to  the  complexity  of 
the  problem.  The  present  work  attempts  to  provide  the  first 
step  towards  the  solution  of  the  problem.  It  is  confined,  how- 
ever, to  the  two-dimensional  case  and  to  constant  temperature 
and  strain  rate  loading  conditions. 

5.2  Model  Representation 

Interlaminar  adhesive  layer  (IAL)  behavior  is  most  conven- 
iently represented  by  the  symmetrical  doubler  model  shown  in 
Figure  16  and  fully  described  in  Reference  16.  Reference  16 
also  provides  the  two-dimensional  linear  elastic  solution  for 
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stress  distribution  of  this  model.  The  corresponding  FEM  net- 
work is  shown  in  Figure  17.  The  analysis  is  focused  on  the 
"boundary"  edge  zone  where  stresses  attain  their  critical  value. 

5 . 3  Plane  Stress  vs.  Plane  Strain  Solutions 

Plane  stress  and  plane  strain  states  may  be  considered  as 
two  bounds  for  the  three-dimensional  solution.  They  represent 
the  situation  close  to  the  free  edges  (y  = ±b)  and  along  the 
midsection  (y  = 0)  respectively.  Stress  distributions  based  on 
assumptions  of  plane  stress  and  plane  strain  for  the  linear 
range  are  shown  in  Figs.  18  and  19.  Comparison  of  the  t and 
o distributions  in  the  two  cases  indicated  only  minor  variations, 
it  may  be  concluded  that  the  two-dimensional  solution  is  not 
affected  significantly  by  the  unknown  stresses  acting  along  the 
third  direction.  The  a distribution  along  y derived  from  the 
plane  strain  solution  (^ig.  20) , (which  may  provide  the  upper 
bound  for  the  actual  three-dimensional  state  of  stress)  attains 
its  maximum  at  the  center  line  (y  = 0) . Generally  it  seems 
that  the  worst  condition  for  failure  would  be  at  the  corners 
(y  = ± b , x = ± c)  . 

5.4  Analytical  Proceudre  at  the  Non-linear  Range 

The  procedure  for  determining  the  stresses  at  the  different 
locations  within  the  adhesive  boundary  zone  under  a given  external 
uniaxial  load  will  follow  the  flow  chart  given  in  Figure  21  and 
includes  the  following  steps:* 

1.  Apply  predetermined  external  a level  at  the  central 
adherend. 

2.  Calculate  effective  strains  and  stresses  of  the  different 
elements  in  the  FEM  network  using  the  linear  FEM  program 
and  assuming  the  same  initial  effective  modulus  for  all 
elements . 

3.  Determine  the  specific  secant  modulus  at  each  element 
based  on  the  experimental  relationship  given  in  Figures 
7 and  8 according  to  the  respective  effective  strain 
calculated  in  step  2. 

4.  Rerun  the  linear  program  with  the  modified  modulus  at 
each  element  according  to  step  3. 

5.  Compare  the  calculated  stress  S^.  at  each  element  as  derived 
in  step  4 with  the  stress  obtained  from  the  experimental 
effective  stress-strain  curve  for  the  respective  calculated 
strain . 

_ 


A full  printout  of  a representative  FEM  program  is  given  in 
Appendix  1. 
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6.  Repeat  steps  3-5  up  to  the  stage  where  the  difference 
between  the  calculated  and  "empirical"  stresses  is  less 
than  2%  of  the  final  stress  value. 

The  different  steps  are  illustrated  in  Fig.  22  for  a representa- 
tive element  k.  The  level  of  convergence  was  achieved  by  sig- 
nificantly fewer  iterations  in  the  case  of  plane  strain  compared 
with  the  plane  stress  solution  (Fig.  23). 

5.5  IAL  Stress  Distribution  at  the  Non-linear  Range 

The  results  were  focused  on  stress  distribution  at  the 
boundary  zone.  The  effect  of  axial  external  stress  cc  (applied 
to  the  central  adherend)  on  the  effective  stress  and  strain  at 
the  critical  point  (located  close  to  the  IAL  edge)  is  shown  in 
Figures  24  and  25.  The  difference  between  linear  and  nonlinear 
solutions  is  pronounced,  especially  in  the  plane  stress  case. 

Lower  effective  stress  and  strain  were  found  in  the  case 
of  plane  strain  compared  with  plane  stress.  The  reverse  trend 
was  found  for  the  shear  (tZx)  and  lateral  normal  (oz)  stresses 
as  a function  of  axial  external  stress  (Figure  26):  namely, 
higher  stresses  in  the  case  of  non-linear  solution  of  the  plane 
strain  state.  In  the  case  of  the  linear  solution,  no  signifi- 
cant difference  is  shown  for  the  two  plane  states. 

Shear  stress  distribution  at  the  boundary  zone  for  oc=53.3  kg/mm2 
is  shown  in  Fig. 27.  Here  again  no  significant  difference  can 
be  distinguished  between  plane  stress  and  plane  strain  solu- 
tions, except  at  the  critical  location  close  to  the  lower  edge 
corner  (see  Fig.  17) . Similar  trends  were  found  for  the  two- 
dimentional  lateral  normal  stress  (Fig.  28) . 

In  most  cases,  the  effective  stress  tends  to  level  off  and 
even  drop  faster  in  the  plane  stress  state. 

The  effective-stress  distribution  at  the  boundary  zone 
(Figures  29,  30)  permit  evaluation  of  the  ductile  failure  pro- 
cess of  the  IAL.  The  region  of  "viscoplastic  flow"  may  comprise 
the  elements  where  the  effective  stress  or  strain  exceed  their 
prespecified  limit*  as  shown  in  Fig.  8.  For  the  specific  load 
level  of  the  present  case,  this  region  is  located  close  to  the 
lower  edge  of  the  IAL  (Fig.  17). 

5.6  Comparison  with  Simplified  Elasto-Plastic  Solution 

A simplified  elasto-plastic  model  of  the  stress-strain 
relationship  (Figure  7)  has  substantial  advantages  over  the 
more  realistic  non-linear  one,  in  that  it  reduces  the  para- 
meters for  describing  the  complex  inelastic  process  to  two  - 
the  initial  elastic  modulus  E0  , and  the  yield-stress  plateau  S^. 

* 

This  limit  may  be  defined  as  the  level  where  initiation  of 
residual  deformation  was  detected  (Ref.  121]). 
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It  may  also  permit  inclusion  of  the  strain-rate  and  temperature- 
dependence  in  the  analysis,  provided  their  effect  on  the  above 
parameters  is  available  [3] . The  respective  stress  distributions 
based  on  this  model  (Figures  31-34)  show  that  while  there  is  no 
significant  deviation  from  the  exact  solution  at  the  critical 
location,  there  is  some  farther  away  from  the  edges.  More  itera- 
tions are  needed  for  convergence  of  the  elasto-plastic  solution 
compared  with  the  more  exact  non-linear  one  (Figure  35) . 


6.  SUMMARY  AND  CONCLUSIONS 

The  present  report  comprises  two  main  parts,  the  first  of 
which  deals  with  bulk  and  in-situ  characterization  of  an  adhesive 
material  and  provides  the  empirical  basis  for  the  second  part, 
devoted  to  analytical  solution  of  the  stress  distribution  in 
the  IAL  at  the  non-linear  range.  The  following  conclusions  may 
be  drawn : 

(a)  The  stress-strain  relationship  of  the  bulk  adhesive 
under  a combined  state  of  stress  may  be  derived  from  the  stress- 
strain  relationship  obtained  under  uniaxial  loading.  This 
follows  the  "effective"  stress-strain  approach  and  von  Mises' 
criteria  for  inelastic  behavior. 

(b)  The  shear  stress-strain  relationship  of  the  adhesive 
in-situ,  its  Young's  modulus  and  ultimate  strength,  may  be  roughly 
considered  as  invariant  for  thicknesses  about  0.2mm;  in  these 
circumstances,  bulk  data  provide  the  basic  parameters  for  a pre- 
liminary assessment  of  the  relevant  mechanical  behavior  of  the 
adhesive  within  a bonded  structure. 

(c)  The  finite-element  method  was  found  to  be  adequate 
for  determining  the  stress  distribution  of  the  IAL  at  the  non- 
linear range.  At  high  external  loading  levels,  non-linear  be- 
havior was  found  to  predominate  in  narrow  boundary  zones  close 
to  the  IAL  edge,  whereas  most  of  the  IAL  remains  at  the  elastic 
(linear)  range. 

(d)  The  non-linear  FEM  solution  serves  for  assessment  of 
IAL  failure.  The  latter  is  manifested  by  initiation  of  plastic 
flow  of  the  adhesive,  which  can  be  related  to  a specific  limit 
point  on  the  effective  stress-strain  curve.  Under  a given  ex- 
ternal load,  a "visco-plastic"  state  prevails  in  the  boundary 
zone  of  the  IAL  with  stresses  in  all  elements  exceeding  the 
above  limit. 

The  present  study  is  the  first  step  towards  more  general 
non-linear  analysis  and  empirical  investigation  of  bonded 
structural  systems.  Future  research  will  involve  time-dependent 
non-linear  stress  analysis  of  the  IAL  with  reference  to  the 
viscoelastic-plastic  nature  of  the  polymeric  adhesive;  another 
phase  will  cover  the  influence  of  environmental  factors  (tempera- 
ture, moisture)  on  the  above  time-dependent  mode  of  behavior. 

_ _ . a 
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GLOSSARY 

half  IAL  width 
constant 
half  IAL  length 
initial  Young's  modulus 
secant  modulus 
tangent  modulus 
effective  strain 
effective  strain  rate 
thickness 

adhesive  thickness 
number  of  iterations 
effective  stress 
temperature 

axial,  transverse  and  lateral  coordinates  respectively 

shear  strain 

shear  strain  rate 

normal  strain 

uniaxial  external  strain 

axial  strain  rate 

Poisson  ratio 

shear  stress 

normal  stress 

uniaxial  external  stress 

nondimensional  axial  coordinate 
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IAL  - interlaminar  adhesive  layer 

CAL  - central  adherend  layer 

EAL  - external  adherend  layer 

FEM  - finite  element  method 


TENSILE  STRESS  CT  (kg/mm2) 


TENSILE  STRAIN  K (•/.) 

FIG.  1 TYPICAL  STRESS-STRAIN  CURVE  AND  THE  RELATED  POISSON  RATIO 
FOR  BULK  EPOXY  RESIN  UNDER  UNIAXIAL  TENSION. 


POISSON  RATIO 


TENSILE  STRAIN 

TENSILE  STRESS-STRAIN  CURVES  OF  EPOXY  RESIN  AT  DIFFERENT 
STRAIN  RATES. 
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STRESS  AND  INITIAL  YOUNG'S  MODULUS  OF  BULK  EPOXY  RESIN 


PHOTOGRAPH  OF  TUBULAR  EPOXY 
SPECIMEN  WITH  SPECIAL  DEVICE 
FOR  MEASURING  ANGULAR  DISPLACE 
MENT  UNDER  TORSIONAL  LOADING. 
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THE  EFFECT  OF  SPECIMEN  LENGTH  ON  TORSIONAL  DEFORMATION  - 
PROCEDURE  FOR  CORRECTING  SHEAR  STRESS-STRAIN  RELATIONSHIP. 


AVERAGE  SHEAR  STRESS  r (kg/mm2) 
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t=1mm 

.•correction  jf  S corrected  curve 
; curve  (bulk)  (bulK) 


/ 
l / 


' 

/ / 


experimental  curve 
(bulk) 


3 ■ -1 


7=  69*10  min 
T = 23°C 


SHEAR  STRAIN  7 (7.) 


FIG.  6 PROCEDURE  OF  DERIVING  "TRUE"  SHEAR  STRESS-STRAIN 
RELATIONSHIP  BASED  ON  TORSIONAL  TEST  OF  TUBULAR 
SPECIMENS  OF  DIFFERENT  GAUGE  LENGTHS. 


EFFECTIVE  STRESS-STRAIN  RELATIONSHIP  AND  RELATED 
PROPERTIES  FOR  EPOXY  RESIN  USED  AS  ADHESIVE  LAYER 


EFFECTIVE  STRESS 


prespecified 

viscoptastic 

limit 


e = 0.0385  min' 


1 EFFECTIVE  STRAIN  e (7. 

THE  EFFECT  OF  STRAIN  ON  THE  VARIATION  OF  POISSON'S  RATIO 
TANGENT  AND  SECANT  MODULI  OF  EPOXY  RESIN  UFED  FOR 
ADHESIVE  LAYER. 


POISSON  RATIO 


TORSIONAL  SHEAR.  STRAIN  (•/.)  ‘ 


FIG.  9 TYPICAL  SHEAR  STRESS -STRAIN  CURVES  FOR  THE  ADHESIVE  LAYER 
IN-SITU  COMPARED  WITH  ITS  BULK  EPOXY  REFERENCE  (FIG.  6) . 


AVERAGE  ADHESIVE  SHEAR  STRESS 
AT  FAILURE  (kg/mm2) 


rf 
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PHOTOGRAPH  OF  DEVICE  FOR  COMBINED  LOADING  OF  ADHESIVE 
LAYER  IN-SITU 


TYPICAL  SHEAR  STRESS-STRAIN  CURVES  OF  ADHESIVE  LAYER 
IN-SITU  UNDER  COMBINED  SHEAR  AND  AXIAL  NORMAL  LOADING. 
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FAILURE  ENVELOPE  FOR  ADHESIVE  LAYER  IN-SITU  (LOADED 
UNDER  COMBINED  SHEAR  AND  AXIAL  NORMAL  STRESS) . 
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SHEAR  STRAIN  y (7.) 

FIG.'  15  THE  EFFECT  OF  HYGROTHERMAL  HISTORY  ON  ADHESIVE 
STRESS-STRAIN  BEHAVIOR. 


FIG.  17  ILLUSTRATION  OF  IAL  NETWORK  AT  THE  BOUNDARY  ZONE 

USED  FOR  FINITE  ELEMENT  NON-LINEAR  STRESS  ANALYSIS. 
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NONDIMENSIONAL  AXIAL  COORDINATE 


AVERAGE  APPLIED  STRAIN  €c  (%>) 


FIG.  25  THE  EFFECT  OF  CENTRAL  LOADING  ON  EFFECTIVE  STRAIN  AT 
THE  CRITICAL  POINTS  IN  THE  BOUNDARY  ZONE  OF  THE  IAL. 


plane  strain 
plane  stress 


AVERAGE  APPLIED  STRESS  CTc  (kg/mm2) 


LATERAL  NORMAL  STRESS  <3o(kg/mm2) 


X=0.9904  X=  0.9942  X=  0.9  981 
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FIG.  28  LATERAL  NORMAL  STRESS  DISTRIBUTION  AT  THE  BOUNDARY 
ZONE  OF  THE  IAL  (NON-LINEAR  RANGE) . 
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(NON-LINEAR  RANGE) . 


EFFECTIVE  STRESS  s (kg/mm2) 


AVERAGE  APPLIED  STRESS  (Tc  (kg/mm2) 

FIG.  31  THE  EFFECT  OF  CENTRAL  LOADING  ON  EFFECTIVE  STRESSES 
AT  THE  CRITICAL  POINT  IN  THE  BOUNDARY  ZONE  OF  THE 
IAL  (PLANE  STRAIN,  SIMPLIFIED  ELASTO-PLASTIC  VS 
I NON-LINEAR  SOLUTIONS) . 


THE  EFFECT  OF  CENTRAL  LOADING  ON  EFFECTIVE  STRAINS 
AT  THE  CRITICAL  POINT  IN  THE  BOUNDARY  ZONE  OF  THE 
IAL  (PLANE  STRAIN,  SIMPLIFIED  ELASTO-PLASTIC  VS 
NON-LINEAR  SOLUTIONS) . 


NONDIMENSIONAL  AXIAL  COORDINATE 


CONVERGENCE  OF  THE  FEM  ITERATION  PROCEDURE  TO  THE 
EXACT  SOLUTION  (NON-LINEAR  VS  SIMPLIFIED  ELASTO- 
PLASTIC  SOLUTIONS) . 
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APPENDIX  1 

PRINT-OUT  OF  TYPICAL  FINITE  ELEMENT  PROGRAM  rOR  STRESS 
DISTRIBUTION  IN  AN  INTERLAMINAR  ADHESIVE  LAYER  WITHIN 
A DOUBLER  MODEL. 


A N 71  ) 


- * - BEST  AVAILABLE  COPY 


CE'’fO  FORTRAN  H 


ILER  OPTIONS  - MA'Tr  F>  a IN  , r PT  = 0?  . I TN  FCf  !T  = 6 C 1 S T?E  = n n 0 . 

souroe  .er'ttc  .noli't  .nodeck ,lc ao.nc"'Ap,ncecit,  in,r. cxprf 
c 

C TUfl  CTPENSIOfl  PLAN  STRATA /STRFSS  f IKTTF  ELEMENT  PROGRA* 

C FCR  GENERAL  f'ATFRTAl 
C 

R r A L * 0 Tin  F ( 9 ) 

C CP MON  NMP.f  EL  .NIMAT  .f.SLC.N  CFT  .rpOOY  .PTYP,  IE  < 200 , 5 > .90  ( 10  ) ,TM  1 0 ) , 
lpT  ( 10  1 .IJT  ( in  ) ,GGT  ( 10)  ,F1  T ( 10  ) .E’T  (1  n ) ,UI2T  ( 1 0 ) ,L?1T  ( in  ) ,P1?T  (1  m . 

?Cm.T).O''(T,?,10).E.U.GG.El,F2.Ul?tL?i,GlP.GOT(3,3,ip), 

Tx(?oo),rono),uLX(Rnn),vLT(?nr) , kooe <200),Tcr(70),^sc(20), 

H S L R T R V (20.2)  . c U R T R Y (20.2)  , E F (IP) 

COMMON /ONE /px  (in,int,G(LO).p(’.10)«C<’.3)»PT(’.f-)>XO(5).VG(5) 
COMMON  /TWO/  I BANC  « NFO  * RI400),  AK(4C0,50> 

CATA  w A TEL  • WAXNP,  "flxPAT,  "AXnW,  vflXSLC 
l > ?on,  ?nn,  10.  50.  20/ 

999”  READ  lOO.NPPOF.  ( T T T(  E ( T ) . T = 1 . <3  ) 

IF  (A'PRQR.LF  .0)  GO  TO  RES 
in2n  PRINT  700  .f-IPRCB  . (TITLE  ( I ) . 1 = 1 .9  ) 

CALL  CAT  A IN  (PAXEL.wAXNF»MAXMAT,PAxSLC.TCTrP) 
wflXPCF  = 2 * P A X N P 
P AX  P I F = o 
DO  1 1 = 1 .NEL 
DC  1 0 = 1. M 
DC  1 K = l.*t 

LL=  T ARS ( ir { j , J ) - IE(i,k)) 

IF  ( LL  .GT  .1*  AXDIF  ) PAXPTF  = LL 
1 CONTINUE 

TRANE  = ?♦ (PAXPTF  + 1 ) 

NEC  = ? * N N p 

IF  ( IRANC.pt  . P A X n L ) GO  TC  c00 
IF ( TSTOP .0 T . 0 ) Gr  TO  °9” 
call  ASEPFL t TSTOF ) 

IF  ( ISTOR.GT.01  GO  TO  = 99 

CALI-  PANSGL(l.AK,R.NEO.IDAfn.MAXCCF.NAXPW) 

CALL  p A Nc  C L ( 2 . A k . R . N'  E C . IRAnO.'AAXECF  . p A x p v,  ) 
frtnt  yno,  ( t , r < 2* t-i i ,r i 2* i » , 1 = 1 , nnf ) 
call  stress 

GO  TO  QROC 

son  PRINT  901.  TPAND.  PAXRW 
GO  TO  9UR9 

ion  FCRPAT(  T5.’x.RAfl  ) 

?C0  FORMAT  ( /PHI  PROBLEM , Tc , AH. . .eA0  ) 

ACO  FORMAT ( 37H1PHTP!!T  TABLE  1..  NCCAL  DISPLACEMENTS  // 

1 13X.4HN0PE.  RX.  1 1 *-  L = X -O I SP  . , gx  , 1 1 H V r Y-riSP,/ 

2 (5x.T12.2E20. 9)  ) 

9 C 1 format ( ///l on  PANEL irTL  = » I 4 » 2 5H  EXCEEDS  VAX,  ALLCLARLF  r,I4// 

1 A OH  GO  ON  TO  NEXT  PROBLEM  > 

geo  CTCP 
Ef  C 


none i onn 
000C2000 
nooCAOPO 
ooncRono 
nopcfono 

OOOCFOOO 

nonooop 
none  7500 
nooepooo 
noncsono 
onni coon 
pool  loon 

P001700P 

00013000 

poni 7000 
nnoipooo 
oopisoon 

0 0 02  C 0 0 0 
00021000 
"0024000 
00025000 
CCO2GGC0 
00027000 
OOP2FOOO 

nno? c ooo 

000’ COOP 
000’ i ooo 

00032000 

non’ 3000 

POO’4000 
COO’5  OOO 

roo’f ono 
noo’soon 

00092000 

00093000 

00095  0m 

00097000 
C0O4F0OP 
00099000 
00050000 
00051000 
00052000 
0005 ’ OOO 
”0059000 
oop'5000 
POP5EOOP 
00057000 
onpePGO0 
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BEST  AVAILA"' c 


Ui.L 


copy 


AN  7’  ) 


cE'Tgo  fortran  h 


HEP  CPTTCflS  - NAI«F  = I"  A T A , CPT=OP  . L IN  ECMT=60  . S t ZE  = 0 0 0 nK  , 

soiiRrr  .rprrrr.NriLTPT  ,n cheek , loao .nop ap , noec it  .tc.ncxpff 

ECPRCI.T  TNE  **  f TFRF  OOOCIOOP 

COPPON  MNP  « ATI  , MN  A T , N SI.  F , N C F T , NROOY  . FT  YP  . T F ( ? 0 n . 5 ) . RO  ( 1 0 ) . T V (1  0 ) . 00002000 
iET(io),i)T(in),RGTn'H,FTTnc>.E?T(n).ui2Tnn),L'?jT(im.n?7Mn),  00003000 
2CC  <3, 1 ) ,QP<.*,»  ,10)  .F.U.CG,  El  .F?,(J12,t.21  .G12.C0T  (3.?,  10  | , COOC4GOO 

3y(200),Y(P0P) , I.ILX ( 200 ) .VLY (20C  I ,KOnr (200)  ,TCC(70)  ><.SC(20)  , OOOCcOOO 

USt'RTRV  (20,2)  . FIIRTRY  (20 .2)  ,EF  * 1 0)  00006000 

CO  250  TMAT=1,N*AT  OOOCOOOO 

f=et(T'iat)  ooocpooo 

L=LT(IMAT)  OOOCOOOO 

GGrGGT(TMAT)  000)0000 

E1=E1T(IMAT)  00011000 

E2=E2T ( IMAT)  00012C00 

L12  = U1PT(IMAT)  0001  2 000 

L 2 1 =U2  1 T ( I P AT ) 00014000 

G12=G12T ( IPAT ) 00015000 


CO 

350 

1=1  . 

* 

00016000 

CC 

350 

J = 1 . 

3 

00017000 

750 

CG ( I . J ) 

= QQT( I , J 

. IP  A T ) 

oooieooo 

IF 

( E 

.EQ. 

0. 

.ANC.GG  .EG.  0.  ) GO  TO  500 

00015000 

IE 

( E 

.EQ. 

0. 

) GO  TO  300 

00020000 

IF 

( l 

.EG. 

0. 

) GC  TC  200 

00021000 

IF 

(GG 

.EQ. 

0. 

) GG  = E*0.5/ (l.-HJ) 

00022000 

GC 

TO 

400 

00023000 

200 

IF 

(GG 

.NE. 

0. 

) U = 1.  - E*0.5/GG 

00024000 

CC 

TO 

400 

00025000 

;on 

IF 

( U 

.EG. 

0. 

) GC  TO  400 

00026000 

1200 
25  n 


F = 2 . * ( 1 . +U ) « GO 


<400 


050 


5C0 


700 

OCR 


IF  (rtOPT.EC.1  ICO  TC  450 
C<3  ( 1. 1 ) = If  ( 1 ,-U*L  ) 

CCI2.1)  = U*OG (1.1) 

CC ( 3 , 7 ) = GG 
G C ( 2 « 2 ) = 00(1  .1  ) 

GC  TC  1000 

CC  ( 1 .1  >=E*  (1-1  )/(  (1+U  )*(  1-2*11)  ) 
CO  (2. 1 )=E*'J/  ( ( 1-HJ)  * < 1-2»U)  ) 

GC  ( 2. 2 ) =0(3  ( 1 .1  ) 

CG ( 3 , 2 ) =GG 
GC  TO  1000 
IF  ( El  • EQ . f. 

IF  ( L 1 2 , FC . 0. 

IF  ( 121  .EC.  0, 

GO  TO  000 
L 12  - U21*E1/F2 


, AN  0 . F 2 . E 0 . 0. 

) GC  TC  700 
) U21  - I 12«F2/E1 


) GC  TC  1000 


CU1.1I  = El/ll. -L12»U21  ) 
GC ( p. 1 ) = U21 » CO (1.1) 
GGI2.2)  r F2/(1.-L17*U21 ) 
GC  ( 3,  ’ ) = G12 
1000  CC  1200  1=1.3 
CC  1200  j=1 , 1 
C('(I.v,.TdAT)=CO(I,J) 

G ('(.), T,IMAT)=GG(I.v) 
continue 
CCNTINur. 

RETURN 
ETC 


00027000 

O0O2P000 

00025000 

0003  C 000 
00031000 
00032000 
00033000 
00034000 
00035000 
00036000 
00037000 
00036000 
0003c000 

0004  0000 
0004 1000 
00042000 
00043000 
00044000 
00045000 
00046000 
00047000 
0004P000 
o 004  5 000 
O0O5COO0 
00051000 
00052000 

0005  3000 
00054000 
00055000 
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JAN  7?  ) CS'360  FORTRAN  F 

FILER  CPTICNS  - MAVE  = V A IA  , CFT  = 02  , L TA  ECAfT  = 6 G , S I2E  = opOOK  , 

soup  r f iPrcctf , r oltft .aooeck  .lcac.ncv  ap,  Loren .ic.acxpff 

SLBROLT  TNE  OAT  A TN  ( « A X E I . V AXfIP  . V A X v AT  , » H ''St  C . I ST  CP  ) 

COP  "ION  MNP  »AFL  .NMAT .A  SL  C.A  CFT .APODY , VJ yp , IF (20n  »5 ) ,RO ( 1 0 > .TP ( 1 0 ) 
1ET  < 10)  ,'JT  { 1 0 ) ,GGT  <1  o ) .FIT  ( 10  ) ,E2T  < 1 0 ) .U12T  ( 1 0)  ,L?1  T ( lc  ) .C-12T  { lr  ) 
2C0 (3.3  J .Q^I 3,3 ,10) .E.U.GG.El .E?,U1?.LP1.G12,CQT( 3.3.10) . 

3X (200) , Y (200) .ULXI200 ) ,VLY (por ) .HOPE (200). I«CC ’0). JSC (20). 
4SURTRX (20,2) , SUPTRY (20.7) .EF (1 0) 

I STO°  = 0 

READ  1 .MNP. PEL , MM AT , A. SIC .A CFT . NPOD Y 

RP IMT  1 00 . AINC .NFL.AmAT , AJSLC.MCPT .MPOCY 
IF (MNP. I F.VAXAF ) GO  T f 201 
ISTOF  = ISTCF  + 1 
PRINT  251,  VAXAJP 

201  IF(MEL .LE.MAXEL ) GO  TC  202 
ISTOF  = ISTCF  + 1 

PRINT  252.  VOXEL 

202  IF (NVAT.LF .VAXVAT > GO  TC  203 
TSTOF  = ISTCF  + 1 

PRINT  253,  V A X v AT 

203  IF (NSLC.I E .MflXSLC ) 60  TC  204 
ISTOF  = ISTCF  + ) 

PRINT  254,  VAX  SLC 

204  IF(ISTOP.EO.O)  GC  TO  205 
PRINT  255.  I ST  OP 
STOP 

205  REAP  2,  (RO( I ) ,TH( It , 1 = 1 .A  VAT ) 

PRINT  101 

PRINT  5l»  ( I ,PO( I ) ,TM I ) , T=1 ,AVAT ) 

READ  10,  (ET(  I ) .UT<  I ) .GCT  ( T ) ,E1  T(  I ) .E2T  ( I ) ,L'12T(  I ) .L21T  ( I ) . 
1G12T(T).I=1,NVAT) 

PRINT  SO,  (I ,ET( I) ,NT( I ) ,GGT (I ) ,E1T( I ) ,E2T( T ) ,1  12T ( I ) ,U21T ( I ) , 
1G12T(I).I=1,N»AT) 

CC  150  IVAT=1.NVAT 
ro  150  1=1,3 

READ  li,  (GOT ( I « J « I"AT) , J=1 « 3 ) 

PRINT  Si.  (OCT ( I , J , IVAT ) .„  = J . 3 ) 

150  CONTINUE 

CALL  VftTERP 
CC  550  IVAT=1,NVAT 
CO  550  1=1,3 

PRINT  82.  (GV( I , J, I«AT ) , J=1 ,7) 

550  CCATINIIE 
FPINT  103 

N=1 

5 READ  3,  V.KOrF (V) ,Y(V ) . Y( V ) ,ULX(V) ,VLY(V) 

IF  ( W-A! ) 4 « E ,7 
4 PRINT  105,  V 

PRINT  52,V,t<OrE(V)  , X (f  ) , Y ( V ) ,L'L  X ( v ) , VLY  (R  ) 

I ST  OF  = ISTOF  -tl 
GC  TO  5 

7 OF  = V + 1 - f 

RX= (X (V)-X (A-l» ) /CF 
RY= ( Y (M) -Y  < N - 1 ) ) /CF 

8 K0OF ( N ) =0 


onoci ooo 
. noocpooo 
, 0OOC300C 
OOOC  3500 
O0OC40O0 
O00C5000 
noncEonn 
O00C700O 
OOOCSOOO 
OOOCSOOO 
00010000 
0001 3000 
00014000 
00015000 

rooiE ooo 

00017000 
OOOlf 000 
OOOISOOO 
00020000 
00021000 
00022000 
00023000 
00024000 
00025000 
0002E OOP 
00027000 
OOOPPOOO 
0002e000 
00030000 
00031000 
00031 010 
00031020 
00031040 
0003 1050 
00031110 
00031120 
00031130 
00031140 
00031150 
0003 1 500 
00031510 
00031 520 
POOM530 
0003 154  0 

nno,,4000 

O00,5000 
0003E  000 
00037000 
O003P000 
n.OTSOCO 
0004C000 
0 0 0 A 1 ooo 

00042000 

O004300C 

00044000 

O0045000 
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BES1  AVAILABLE  COPY 


v ( r ) = x ( n - 1 i +px 

V ( M = Y ( U - 1 ) ♦ R Y 
L'LX  K )=1.0 

VI  Y (A  1=0.0 

6 PRINT  5?«l\  .KPTE  <N)  . V <N  1 . Y (M  .1  LX  <N>  . VLY  (l\  ) 
NrN+l 


IF  (M-iv  >°,6 .« 

9 IFfN.LE.NNPl  PC  TO  *5 


FPINT 

106 

13 

L 

= 0 

14 

READ 

15.  M 

, < TF  ! w, T 1 , 1 = 1 .5  1 

16 

L 

=L*1 

IF (M-L 1117,17,10 

117 

PRINT 

1 1 fi , 

FP  KIT  53.  A'.  <1E(M,I),  1 = 1.5) 

I STOP  = T ST  Cd4 1 

gc  to  m 

ip  ieu  .1 1=  teil-i  ,i  Mi 

IE<L .21=  IF(L-1.?M 
IE ( l , 3 1 = IE ( L - 1 , 3 1 + 1 
IE (L . 4 ) =TE (L-l .4 1+1 
IE  <L ,5)  = TE  <L-1 .5) 

17  FPKJT  53.  L.  < TE<L.  II  . 1 = 1. = 1 
IF ( M-L  120.70.16 

20  IF (NtL-L 1 2 1 .21 . 14 

21  CONTINUE 

IF (NSLC.EG  .0  1 GO  TO  31 

30  fpint  ion 

DC  40  L = 1 .NSLC 

READ  4 1 » ISC  (I  ) .USC(L)  .SURTRXL  ,1)  .SURTRX(L,?1  .FLPTPY(L.l)  , 

1  SLRTR Y ( L » 2 1 

40  PRINT  4? . TSC ( L 1 , J5C (L 1 « SUP  T R Y (L.l ) . SLRTRXIL.21  , SUPTRY (L . 1 1. 
1 SURTPY(L.P) 

3 1 IF(ISTOP.EO.O)  GC  TO  555 

PRINT  900.  ISTOP 


1 FORMAT ( 615 1 

100  FORMAT ( 35R0INFUT  TAOlE  1..  pA'IC  PARAMETER?  // 

1 5X.  4 01-  MUMPER  OF  A C T A L POINTS. Tr  / 

2 5X  • 4 01-  MUMPER  OF  FLFvFrTS 

3 5 X . 4 p P MUMPER  OF  rTFcpPENT  MATERIAL' ,1  = / 

4 5X,  4 OP  MtjMREp  OF  5LRFACE  LOAD  C-ARCS , 1 = / 

5 ax,  4 0 P 1 = PLATE  STRAIN,  2 = PLANT  'TRESS.  . .,1  = / 

6 *x,  4 0 P POCY'  FORCES ( 1 = IN  -Y  OIREC..  0 = NONE) , T 5 1 

251  FOR  "1  AT  (////■*  ?p  Tor  '*  A P Y NOCAI.  POINTS,  MAXIMUM  = , 15) 

252  FORMAT  (////’OP  T^C  "Af.Y  EL  rMrP  TS  , MAXK'OM  = ,1=) 

25  7 FOR  NAT  ( / ///30P  TOC  "'ANY  MATERIALS.  M A X I vl  M ■=  , T c 1 

254  FORMAT  ( ////40P  TOC  "ARY  ' L R F -A  C E LOAr  CARES,  MAVKHM  = ,1=) 


25  = FORMAT ( ////2PP  FXFCUTICN  pAl7Fr  BECALSE  CE.T5.l3M  F A T A | rFRCFS/1 
2 FORMAT ( 2E1 0 . 3 1 

1 C 1 FORMAT ( 36HnifFUT  T " nL  E 2..  MATERIAL  PROPERTIES  // 

1 10P  wfTERIA|  , c X , 1 0 P v 9C UL MS  OF , 6X . 9HP0 I 'SON • S . 7X « 

2PP  MATERIAL  OPMATEPIOt  / 

34X  ,FiHNUmBFR«5X  .1  OPEI  ASTTCTTT.0V.7M  RATTC.PX«7(-Crf  SITY.RV. 
45PTMICKNESS  ) 

5 1 FORMAT ( I10.2E15.4) 

10  format { ml  1 0 . 7 ) 

0 0 FORMAT  ( T 1 0 , or  1 5 ,4 ) 


00046000 

00047000 

00046000 

00045000 

00050000 

non1 2 3 4 5 6 i ooo 

00052000 
00053000 
00056000 
00057000 
000=6000 
0C055C00 
00060000 
00061000 
00062000 
00063000 
00064000 
00065000 
00066000 
00067000 
00066000 
00065000 
00070000 
00071000 
00072000 
00077000 
00076000 
00077000 
00076000 
00075000 
000=0000 
000=1000 
00062000 
00063000 
00064000 
00065000 
0006  6 0C0 
00067000 
00066000 
00065000 
00050000 
00051 000 
000=  2000 
00053000 
000=4000 
0005=000 
000=6000 
OQO570PC 
00056000 
000==000 
nm  crooo 
noi r i ooo 
001C2000 
00103000 
001C400C 
001C5000 

nnir=oio 

001C=  020 


60 


BEST  AVAILABLE  COPY 


1 1 FORMAT  ( 3E1  fl . 3 ) 

61  FORMAT (3 Fin. 3) 

6?  FORMAT (3E10.3) 

103  FORMAT  ( 34H1  INPUT  TABLE  3..  NOCAL  POINT  0 AT  A // 

1 5X,5HNCnAL«4RX.7HX-CI':F  . ,PX,7l-Y-DISr./ 

25X.5Hn0lNT,6X,4HTYPE.14X,l(-X«l4X*lHY,ex.7H0R  LOAD.  PX  , 7 HOP  l f AT  ) 
3 FORMAT (2I5«4E10.3) 

1CC  FORMAT (?X,17HERR0R  IN  CART  NO. ,15/) 

3?  FORMAT (2110. 4F15.4) 

10  6 FORMAT  ( 34H1  Tr!FUT  TABLE  4..  ELFMEMT  DATA  // 

1 1 IX  , 316 GLOBAL  IN  DICES  OF  ELEMENT  NODES/’ X , 7RELEMFNT , 

27X.1H1 ,7X,lH?,7X,lt-3,7X,lF4,0>,pHIMATERIAL) 

IIP  FORMAT!  5X.  25  HFRRCR  IN  ELEMENT  CAPO  NO.,Ic  ) 

15  F ORMAT (615) 

53  FORMAT (T10 .4TP , 110 ) 

10P  FORMAT ( 37H1INFUT  TAPLE  5..  SURFACE  LOADING  DATA  // 

1 1 7 X « 3 3HSURE  ACE  LOAD  INTENSITIES  AT  NODES/ 

24X  , 6HN00E  I « 4X  , PHN’CnE  J « 1 OX  , 5R  X I » IPX  . 2HXJ  » 1 OX  , ?f- Y I . 1 OX  , 2hY ) 

41  TORMA,  (2I5.4E10.3) 

42  FORMAT (2110 .4E 12.4  ) 

ROD  FORMAT  ( ///45H  ASSENPLY  Al\ir  SOLUTION1  LlLt.  NOT  BF  PE  R F 0”  met  * , I 5 . 

1216  RATAL  CAFD  ERRORS  ) 

999  RETURN 
EN-C 


onic'030 
001 C' 040 
001C5050 
O01CF000 
001 C7000 
O01 CPOOO 
001C9000 
OOllCOOO 
oon  iooo 
00112000 
00113000 
00114000 
00115000 

o 0 1 1 £ 0 0 o 

00117000 
P011P000 
00119000 
001 20000 
00121 000 
00122000 
00123000 
00124000 
00125000 
00126000 


o A N 73  ) 


C S / ’ 6 0 FORTRAN  H 


FILER  OPTION'S  - NAMFr  A T N , CFTr  0?  , L IN  ECNJT  = 60  . S I2F  = o 0 0 OK  , 

SOURCE  tEBCD  TC  .NOLISI  .NDDECK  ,LC  AD  .MCM'AP  , NOECIT  . ID  iNCXPEF 
SUPROLTINJE  orc«PC(L,N) 

COMMON.  / TWO/  IF  AUO  < N EC  .F  ( 400  ) . AK(400*C0> 

C THIS  SUBROUTINE  MODIETFS  THE  A'REwPLAGE  STIFFNESS  ANir  LCAOP  FCR  ThF 


FRESCRIBED  DISPLACEMENT  L AT  DEGREE  OF  FREECCM  N.  EQ.16-19B).  (REF, 
DC  100  Mrj.lRAND 

K = N - M + 1 
IF ( K . LE . 0 ) CO  T 0 50 
R ( K ) = R • M ' - AK(K,w)*tj 

AK ( K , M ) =0.0 
50  K=N+«-l 

IE(K.OT.NFO)  CO  TC  ICO 
R ( K ) = R ( K ) - AK(N,W)*U 


AK(M.M)  = 
ICO  CCNjTINUr 

AK(M.l)  = 
R (N  ) = U 

RETURN 

END 


0.0 


1.0 


OOOC 1 000 
000C200Q 
C0OC300O 
l)O00C4000 
O00C500C 
00006000 
O0OC7O0O 
OOOCPOOO 
OOOCSOOO 

oooirooo 
0001 1 000 
00012000 
00013000 
00014000 
00015000 
00016000 
00017000 

nooi pooo 


k L 


> 


AN  72  ) CS  /?60  FORTRAN  P 


I LER  OPTIONS  - NAPPE  = P'1  A TN  , OPT  = 0?  . L IN  ECNIT  = 6 0 , S 1 2E  = A 0 0 OK  , 

SOURCF  , FRCD  TC  « P PL  T S T « N OOECK  « L C AD  , NOP*  AP  , N OFC  I T , I D » N C X RF  F 


ccp'mon 

NNP  .NFL  .PIMAT  .NSLC.P  CFT  , piPODY  , 

P'TYP.IF(20(',5).R0<10),TM10), 

00012000 

1FT (10) 

,UT(10),GGT(10),ElT(10).E2T(10),Ul2T(10),L?lT(10),rl2T(in), 

00OC2000 

2CC<3.3 

).op»(3,;.i0).r.u.GG.ri«F2,ui2, 

L21 .GlP.QOf (3.2.10), 

P10C2500 

3x ( 200  ) 

, T ( 2 0 0 ) . IILX (2  00  ),VLY(20r),KOrE(2QO),I  = C<20),v.SC(2P), 

OOOC400P 

4$L'RTRX  (20,2)  .SURTRY  (20.2)  ,EF  ( 1 0) 

00005000 

COP' MON /ONE/  ox  (10. 10),  0(10).  ”(3.  10), 

C(3*2).RT(?,6),X0(5).YG(f) 

00DC600C 

CCP'MON/TWO/  TEAND,NEG.R(40C) . AK(400 

.50  ) 

OOOC7POO 

riP'ENSlON  LP  ( 6 ) 

OOOC5COO 

nooCFooo 

REWIND 

1 

oooi coon 

TSTPF  = 0 

00012000 

RT ( 1 ,4 ) = 0.0 

00014000 

RT(1 ,5)  = 0.0 

00015000 

PT ( 1 ,6  ) = 0.0 

00016000 

RTI2.1)  = 0.0 

00017000 

RT  ( 2 . 2 ) = 0.0 

00015000 

PT  ( 2 . 3 ) = 0.0 

00015000 

DC 

2 1=1. NEQ 

00026000 

R ( I > = 0.0 

00027000 

DC 

2 J =1 . IRAN  0 

00026000 

2 

AK(I.J)  =0.0 

00025000 

DO 

10  M=1 ,NEL 

00022000 

IFlinw.N)  .RT.O)  GO  TO  11 

00024000 

ISTOF  = TSTCF  + 1 

00025000 

GO 

TO  10 

00026000 

11  call  QIJAP(M,AREA) 

ono’7ooo 

IF ( ARFA.OT  .0.0)  GO  TO  16 

00025000 

ISTOF  = TSTOP  + 1 

0002c000 

PRINT 

20  ,r 

O0C4  0000 

16  IF(IE(P».3>  .FQ.TF(P',4)  ) GO  TO  26 

00044000 

DC 

31  J = 1.2 

0004=000 

IJ=  10-J 

00046000 

TK=  1J+1 

00047000 

PIVOT  = OK ( IK  . IK  ) 

00045000 

DO 

32  K = 1 . I J 

0004=000 

F = GK ( IK ,K ) /PIVOT 

00050000 

’ GK ( IK , k ) =F 

000=1000 

CC 

33  T = K , I J 

00052000 

OK ( T ,K )=OK ( T ,K  ) - F»OK(T.IK) 

000=3000 

2 3 

OK ( K , T ) = OK ( I ,K  ) 

000=4000 

2? 

0 ( K ) =0(K ) -GK ( TK ,K ) *C ( TK  ) 

000=5000 

21 

0 ( I K ) =0(IK)/PIVCT 

00056000 

26  WRITE 

(1)  ((OK(I,J),J=l,10),T=q,10), 

C ( R ) . 0(10). 

00061000 

1 ( (0 ( I . 

J),J=1.10),T=1 ,3), ( (C(I,J)«J=1 

« 3 ) . 1 = 1 .3  > . XQ ( 5 ) . YC(  = ) 

0O062000 

L I P'  = P 

00066000 

IF  ( TF  ( m , 3 ) ,F0,IE(P'.4)  ) LT"  = 6 

00067000 

00 

40  1 = 2. LIP". 2 

00066000 

IJ  = 1/2 

0006=000 

LP(I-l)  = 7*TF(P’.T^)  - 1 

0O07C000 

40 

LP(T)  = 2*TF  (P'.IJ) 

00071 OOC 

DC 

50  LL  = 1 .LIP' 

00072000 

I = IP'LL ) 

00073000 

R ( I ) = R ( I ) + 0 ( L L ) 

00074000 

DC 

50  P"»=1  , L IM 

0007=000 

J = IP(P'P')  - I 1 

00076000 

62 


BEST  AVAILABLE  COPY 


IP  'J.LF.  0 ) GO  TC  50 

AK(T.J)=  AK(T,J)+  GX(IL«MI“) 

50  CC’v'TIMIF 

10  CCNUMIJP 

DC  65  N=1«NNP 

if (kocf (Nv.ro.3 » go  tt  cc 

Kr?*N 

IP ( KCnr (A ) .EG .1 ) GO  TC  67 

R(K-1  ) = R(K-1 ) + I L* ( N » 

IF (KOCE (M) .NE.O ) GO  TO  5C 
5 7 R <K ) = R <K  > + VLT (A  ) 

5*  CCFlTINUF 

IF  (FiSLC.FG  . 0 ) GO  TO  60 
00  61  L = l.NSLC 
I = 1 RC  ( l ) 

J = wSC(L  ) 

T T = ?«  I 
*J  J “ 2 * 0 

OX  = X ( J)  - X ( T ) 

CY  = Y ( J ) - Y t T ) 

FL  = SORT(OX*nX  + CY»PY) 

FXT  = SURTRX (L  . 1 ) *FL 
PXJ=SURTRX (L.2)*FL 
PY I=5URTRY (L . 1 )*F L 
PYJ  = CJRTRY (L .?)  *FL 

R(TI-l)  = R(IT-l)  * FXT/3.0  + PXJ/6.0 
R(JJ-l)  = R { J J - 1 ) ♦ F X T /6 . 0 + PXJ/G.O 
R(TI)=R(TI)4  PYI/2.0  + PYJ/6.0 
R(JJ)=  R(LJ)  + PY 1/6 , 0 + PYv/^.C 
€1  CONTINUE 

60  OC  TO  M=J ,MMp 

IF  (KOCEO“>  .GE.O.AM1.KCCF  (I*  ) .LF.3  ) GO  TC  72 
ISTCF  = ISTCF  + 1 
GO  TO  70 

7?  IF(KOOE<w> .FQ.O)  GO  TC  70 

IF(K0pF(|w).FQ.2)  GO  TC  71 

CALL  GEC^FC  (IJLX  <R  ) > 

IF (KOOF (M)  .EO.l ) GO  TC  70 

71  CALL  GECRFC  < VLY  (F  ) ,?*M 

70  COffTIMir 

FNCFILE  1 

IF ( TRTOP .FO.  0 ) GC  TO  FI 
PRIMT  100,  IGTOP 

20  FORMAT  (/5X.  1 ■’F  ARFA  CF  FLFRFA'T  ,15. IRE-  I o 1 F G A T I V F /) 

100  FORMAT (////R?F  SOL  L T I0f  WTLl  ACT  PE  FERFGRYFn  npcALSE  OF  , I * , 
1 15H  CATA  fpponc  /) 

01  RETURN 
FNC 


00077000 

P007POOO 

nousooo 
no  of coor 

OOOf ’OOP 
OOOfMOOO 

nooefooo 

ooneeooo 

O00P7000 

ooof eooo 

ooopfooo 

00050000 
0 0 06  <4  00  0 
00005000 
00006000 
00007000 
OOOOPOOO 
ooooooon 
001 coooo 
omciooo 
noicpooo 

001C7000 

OClCROOO 

OOlCfOOO 

oniceooo 

001C7000 
001 CPOOO 
ooicoono 
non  ooon 
coniooo 
00116000 
00117000 

oon  eooo 
00116000 
00120000 
00121000 
00152000 
POl2*000 
001*4000 

noi2=ooo 

noisgooo 

00127000 
OOIGPOOO 
''O12F000 
OOlJCOOO 
pot » 1 000 
P01Y2000 
noiY?ooo 


63 


BEST  AVAILABLE  COPY 


'M  7 3 


C5'360  fortran  H 


LER  OPtTOMS  - NIAMEr  MAIN ,CFT=02 .L TNECUT=60 , SIZE=OOOOK , 

RCUPPE  , FBCC IP  ,MOl.TST  .NCDECK. LOAD. NOMAP, MOEC IT  , IC.NCXPFF 
GLRROLTINE  Ol)A  D ( M , TOT  Al  A ) 

C Cl*  MON  hNP  < NFL  . M«*  AT  . N SLC  , N CF  T , NROD  Y . I*T  YP  . IF  (200.5). RO(IO).TMIO). 
1ET ( 1 0 ) , UT ( 1 0 ) « GGT ( 1 0 ) iFlTdC'  , E2T ( 1 0 ) » U12T ( 1 0)  , U21 T (10)  , Cl  2T  (1(1)  , 
2CG (3,3 ) .O'*  (3,3 ,10) ,F,U«GG,E1 ,E2,U12.L21.G12.QQ7 (3,3,10  ) . 

3X (200) , Y (20A J ,ULX (2no ) , VLY  < 2 0 C ) .HOPE (200  ),T5C(20),uSC(20). 

4SLRTRX  (20.2)  .SURTRY  (20.2)  ,EF  ( 1 0 ) 

COMMON /ONE/  QM10.1O),C<lC>.P(3.10>.C<3.3).PTC'.6),XG<5),YG(‘:) 
CCMMON /T  WO / TEAND.NFG.P(4rQ>  . AK(400,=0) 


1=  TE (M, 1 ) 

J=  IF ( M , 2 ) 

K=  IE  (M.  3 ) 

L=  TE  ( M , 4 ) 

MT YP  = IE ( M , 5 ) 

TOTflia  = 0.0 

IF(fJMAT.EQ.l.AND.w.GT.l)  PC  TC  5 
C(1.1)=PM(1,1 .MTYP) 

C(1«2)=GM(1,2,MTYP> 

C(1*3)=GM(1.3,MTYP> 

C(2.1)=0M(2,1.MTYP) 

C(2,2)=OM(2,?,MTYP) 

C(2i3)=0M(2,3.MTYP) 

C(3.1)=GM(3.1.MTr Pt 
C(3.2)=0M(3.2.MTYP) 

C(3«3)=0M(3,3,mtyP) 

5 LIM  = 4 

IF(K.EG.L)  LIP  = 3 
X Q ( 5 ) = 0.0 
YQ(5)  = 0.0 
DO  TO  N=1  , L IM 

NN  = TE(M.NI) 

X 0 ( M ) = X(NM) 

YQ(N)  = Y ( NN  ) 

XG  ( 5 ) = XQ(5>  + XIAN'I/FLOATILI1*) 
10  YG  ( 5 I = YG  ( f ) + Y (NM  )/FLOAT<(.  IM  > 


DC 

13  II  =1,10 

0 ( t i ) = n . o 

DO 

12  JJ  = 1,10 

12 

QK(TT ,JJ)=0.0 

DC 

13  JJ=1 , 3 

13 

fi(JJ.IT)  = 0.0 

IF(K.NF.L)  go  to  15 


CALL  CST(1,2,3,T0TALA> 

GC  TO  PRO 

15  CALL  C R T ( 1 , 2 , = .AREA) 

Total  a = totala  + area 

CALL  CST(2,*,5  .AREA) 

totala  = totala  + area 

CALL  CRT ( 3 « 4 « = .AREA) 

TOTALA  = TOTALA  + AREA 
CALL  CRT (4,1,'. ARFA  ) 

TOTALA  = TOTALA  4 ARFA 

959  RETURN 
CMC 


onoc i ooo 
00002000 
O00C3000 
00003500 
000C4  000 
00005000 
OOOCE  000 
00007000 
OOOCPOOO 
0000=000 
00010000 
0001 1 000 
00012000 
00013000 
00014000 
00015000 
00020000 
00021000 
00022000 
00023000 
00024000 
00025000 
0002F000 
00027000 
00020000 
00033000 
00034000 
00035000 
00036000 
noo’ 7000 

O0O3P00O 

00035000 

0004000P 
00041000 
00042000 
0004=000 
00046000 
00047000 
00040000 
00045000 
00050000 
0005 1 000 
000=2000 
0C05  3000 
00054000 
"0055000 
000=6000 
00057000 
OP05P000 
0OO==00O 
O006E000 
0006 1000 
00062000 
0006  300U 


I 


64 


BEST  AVAILABLE  COPY 


H -»3  ) 


C'deo  FORTRAA  H 


[ LFR  OPTICUS  - AJAf'Fa  MATA  ,CPT  = n?,l  TAFCNT  = 60 ,ST2F  = 00nnK , 

sonnrF.E:RrrTr,MnLTST.ArnECK,LCAn,A'rvAP.KOFrTT.Tn.ACVRrF 
SURROLTINF  CSTd.o.K.APFA  1 POOPIOOO 

CCMMOA  MNP  ,ATL  , NMAT  .*  SLC .(  C F T , rnoPY , M T YP . T F (20n.3,),RC(10),YMin>,  OOOCPOOP 
1FT  ( 10  > ,UT  ( 10  ) ,GGT  ( 10  ) ,F1  T d 0 ) ,F2T  ( 1 0 ) ,!J12T  <1  0 ) ,L?d  ( 10  ) .C127  (1  0 ) , 000C300P 
?CG  ( 3,  * ) ,0M(  3 .3  . 1 0)  ,F.U  .GG.E1  .E?«IM  ?,L?i  , C-1  2,GGT  ( 3.3 . 10  ) . P00C3S00 

3X(P001,Y(20rn,U|X(200).VLY(?nrt ,KnnF (?nO).ISC('’0).vSC<20).  O00C400P 

4SURTRX (20,? 1 , «UR TRY (20,2) ,EF  d 0 ) OPOCFOOP 

cor*  MCA  /dJF  / OK  (10*10)  ,0(10)  » n ( 3 , 1 0 ) , C ( 3 * 3 ) « RT  ( "*  » 6 1 » X <3  ( 5 ) , YC  C 5 ) OOPC6000 

PPMMOA /TWO/  TFANO.AFC.P <400  ) , AKI4O0.F0)  OQ0C7000 

C IMENS  T ON  CR ( 3 <A) ,LC (6 ) ,LT ( 3 ) ,TK <6,6 > POOCeOOP 

r ooocsooo 

LT (1  I = I 0001  0000 

LT ( 2 ) = J 00011000 

L T ( 3 ) = K 00012000 

PT(1,1)=  YO(J)-YG(K>  00015000 

PTd,2)  = Y 0)  ( K ) - VG  ( I > 00016000 

RT(1,3)  = YGd)-YG(J)  00017000 

PT(2,4)=XG(K)-XQ(J)  O001P000 

eT<2.5)  = XC  (T)-XG(K)  00016000 


PT(2,6)  = YC(J)  -XGU) 

OOC2COOO 

PT(3.1)=PT(2,4) 

00021000 

PT(3,2)  = PT  (2,5) 

00022000 

RT ( 3 , 3 ) = RT (2,6) 

00023000 

RT ( 3 , 4 ) = RT ( 1 , 1 ) 

0 0024000 

RT(3.5)=  RT{  1,?) 

0002=000 

RT ( 3 » 6 ) = F T ( 1 , 3 ) 

00026000 

ARFA  s(ST(9,4)*BTU  .3 

) - PT(2,6)»RT( 1 ,1 ) )/2.0 

0002700 

00 

10  I T = 1 ,3 

0007000L 

00 

10JJ  a 1,6 

00031000 

CR ( T 1 . JJ ) a 0.0 

00032000 

CO 

10  1 

KK  a 1,3 

03033000 

10 

CB(IT.JJ)  = CP(IIu'J) 

+ C ( 1 1 . KK ) *PT (KK , Jo  ) 

00034000 

DC 

12  TT  = 1,6 

00037000 

DC 

12  JJ  = 1.6 

0003P0P0 

TKITl , J J ) = 0 , 0 

000*3000 

OC 

12  KK  a 1 , 3 

00040000 

1? 

TKITl, JJ)a  TK ( T T , Jt 

) +PT ( KK , I r ) *CP (KK , JJ ) 

00041000 

DC 

15  11=1,3 

00046000 

LC(IT)  = 2 *L  T ( T I 1 - 1 

00047000 

1* 

LC  ( Il-t-3)  a 2 *L  T ( T I ) 

00043000 

DC 

30  T T = 1 ,6 

00046000 

LL  = LC(II) 

000=0000 

FK  = 1.0/(4«0*ARFA) 

00051000 

F R = 2,0*FK 

00032000 

DO 

20  JJ  = 1 .6 

0n05?000 

MM  = LC(JJ) 

00034000 

20 

GK ( L L ,MM)  = OK (LL ,M*») 

+ TK(II,Jj)*TP(MTYP)*FK 

0003*000 

00 

30  JJ  = 1,3 

00036000 

30 

P(JJ.LL)  a P(JJ,LL)  ■* 

RT  ( JJ.  I I ) *FF 

O0O57000 

ec  q 


IF (NBCPY .FO.O ) GO  TO  <366 

TRCTYF  = ARF  A * RC(wTYP)* 
ROrYF  = •TPCFYF/3,0 

□ 0 3p  TT  = 1 .3 

JJa  2*  LT < TT  ) 

0(JJ)=  Q(JJ)  + °or  YF 


IF (MTYP) 


PETHRA 

FAC 


ooofoooc 

OCOE1000 
0006  2000 
OOPF  3000 
00064000 
00063000 
00066000 
OOOf 7000 


I I 


BEST  AVAILABLE  COPY 


- 65  - 


AN  7! 


CS'360  FORTRAN  H 


ILER  OPlICMS  - NAFF  = FAIN ,CPT  = 02 .LINECNT  = 6C ,Sl2f  = OOOOK  , 

SOURCF .EBCOIC.NOLTST . N OOECK . LOAD . NCPAP . NOED IT  . IO.NCXREF 
SUBROUTINE  STRESS 

COFFON  NNP  • NFL  . NIF  A T , N SLC  . N CP  T , NpODY  . F T YP  . TF  ( 20  0 . 5 ) . RO  ( 1 0 ) . TK  1 0 1 t 
1ET ( 10) ,UT< 10) , GGT (10) . E 1 T ( 1 0 ) , E2T ( 1 0 > * U1 2 T ( 1 0 ) . U21T ( 1 0 ) . G12 T Cl  0 1 . 
2GQ(3,3).OF(3,!,10).E.U.GG,E1 . E 2 , U12 . L 2 1 . G 1 2 , GOT ( 3 . 3 . 1 0 ) . 

3X (200) «Y (200) « UL  X ( 2 00 ) « VL  Y C 2 0 C ) .KODE1200  1 .ISC<20),JSC(20), 
4SURTRX(20,2l,SIJRTRVf?0.21.EFM0) 

COPFON/ONE/  QF<10.10>.G(101.n(3.10).C<3.3).PTd.6),XQ<5).YC<51 
COPFON/TWO/  IF  AND. NEC  .R (400  ) . AKI400.50) 

C I FENS  T ON  STGI6) 

C 

PEWlNn  1 
PRINT  300 

NOLINE  = 4 7 
DO  5 Prl.NEL 

RE  AO ( 1 ) ( <QK ( I . J) , J = 1 .10)  , 1 = 1 ,2)  , 0(9),  0(10), 

1 ( (8<I« J) tw=l.I0) • 1=1 •?) . <(C(I.J).J=l*3>.I=1.3)t  XC.YC 

LIF  = 4 

IF< TEIF.3) .EQ.IE (P.4)  ) LIF  r 3 
00  10  1 = 1. LIF 
II  = 2*1 
JJ  = 2* IE ( F . I ) 

G ( I T - 1 ) = R< JJ-1  ) 

10  G ( 1 1 ) = R ( JJ 1 

IFUIF.E0.3)  GO  TO  16 
DC  15  K= 1 . 2 

JK  = K + P 
IK  = JK  - 1 
DO  15  L=1.TK 

15  0 ( JK  ) = G ( JK  ) - GK<K,l)*0(L> 

LIF  = 10 

FAC  = 0.25 
GC  TO  17 

16  LIF  = 6 
FAC  = 1.0 

17  DO  20  1=1.3 
EP(T)=0.0 

DC  20  J=l.  LIF 

20  EP<I)=EP(I)+P<I.J)*Q(J)*FAC 

DO  30  1=1.3 

SIG  < I 1 = 0.0 
DO  30  J=1  • 3 

3 0 SIG<  I )=SIG(I  )4CU.  J)*EF  tw  ) 

SP  = I SIG ( 1 1 +SIG (2 1 1/2.0 
SF  r < SIG ( 1 1 -STG ( 2 1 1/2,0 
DS  = SORT (SF*SF+SIG(?  )*SIG(3)  ) 

S I G ( 4 1 = SP  ♦ DS 

s ig  f J i = «p  - ns 

SIGI61  = 0.0 

IF(SIG(3). NE. 0.0. AND. SF.NE. 0.01  SIG<6>  = 20 .64 0 » AT ANp ( S IG  I 3 > . 
1 SF  I 

IF(NOLINE.GT.O)  GC  TO  *4 
PP INT  1 00  0 

NOLINE  = 49 

•4  NOLTNE  = NOLINE  - 1 

5 PRINT  1010.  F, XC. VC. (SIG(I). 1=1,6) 

ENCFILE  1 


000C1000 
00OC2000 
00003000 
00  OC  3500 
000C4000 
COOCSOOO 
OOOCfOOO 
">0007000 

oooceooo 

000C5000 

00010000 

00011000 

00012000 

00016000 

00017000 

oooieooo 
00021000 
00022000 
00023000 
00024000 
00023000 
00026000 
00027000 
00030000 
00031000 
00032000 
00033000 
00034000 
00035000 
O0O3P000 
00039000 
00040000 
00041000 
00042000 
00043000 
00044000 
00045000 
00046000 
00049000 
00050000 
00091000 
00092000 
00034000 
00055000 
00056000 
00057000 
0 009  P 0 00 
00055000 
0006C000 
00061000 
00063000 
00064000 
00063000 
00066000 
00067000 
O 0 06  P 0 00 
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BESTAVAI1ABLE  COPY 


*CO  FORMAT  (47H10I  TPUT  TAPLF  ^TpfSSFS  AT  FLFMFNT  CFNTRCID*  // 

11  X ,7HF|_FMENT  .EX  . lHX,®X.lUY.4X,AHSIGMA<X).4X,PHSlGMA(Y).4y, 

2»HT AU (X . Y ) ,4Y , PHSTGMA ( 1 ) .4X .nPSIGMA ( 2 ) , 7X,*HAPGLF  ) 

10  OO  FORMAT ( 1H1 , 7HELFMFNT,RX,lMX,9X.lHY,4XeHRTGMA(X),4X,PH«iri'A(Y), 
14X  ,0HT  AU<  X , Y ) ,4X  , AI-STGMA  < 1 ) , 4V  , PHSIGM  A ( 2 ) . 7X.PHAP.GLE  ) 

1010  FORMAT (IS.  2F1 0.2.1P6E12.4 > 

RETURN 

ENC 


0006*000 

00070000 

00071000 

00072000 

00077000 

00074000 

00077000 

00076000 


JAP  7J  ) 


CS'760  FORTRAN  P' 


3 I LEP  OPTTCPIS 


PI  A ME  = MAIN  , OPT  = 02  ,l  TNFCP)T=G0  ,SI2E  = OOOOK  , 

R<?l'pCF  .ERCC  IC.f.OLTST  .NOOECK  .LCA0.NCMAF>,NOECrT.  m.NCXREF 


SL'PROLT  INF  0 A N SOL  ( K KK  »AK  .P  « P.  F G . iBAfn.NQlM  « M n I M | 

000C1000 

SYMMETRIC 

PAMC  "ATRIX  EOLATION  SOL  VFR • (REF.  2) 

OCOC?  000 

ooor jooo 

KKK 

r 1 

TRIANGULAR T7ES  TPF  RANG  M A T R T X AK , EG.  (2-5) 

000C40O0 

KKK 

= 2 

SOLVES  FOR  RIGP-T  hand  sidf  r<  solution  peturpjs 

IN  R.  FC.12-3)  000C7300 

POOCf ono 

P I MENS  I ON  AK  (PniM,MnjM) . P C 1 > 

00007000 

PRC  = PIEC  - 1 

OOOCEOOO 

NR  = PIER 

OOOCSOOO 

IF(KKk.fo.2I  GC  TO  200 

0001 0000 

□c 

120  N=  l.NRS 

00011000 

Mr  p>-l 

00012000 

MR  = MImO(TPAPO.PR-M) 

0001 7000 

PIVOT  = AK(N.l) 

00014000 

oc 

120  LrJ.MR 

nooisooo 

CP=  AK  (PJ.L  P/PIVOT 

00016000 

I = M+L 

00017000 

j = r. 

0001 POOO 

DO 

110  k=l ,mr 

0001S000 

J = w + 1 

OOC2CGOO 

11° 

AKIT.JIs  AK(T.J)  -CP«AK(N,K) 

00021000 

120 

AK(M.L)  = CP 

00022000 

GC 

TO  400 

00027000 

?C0 

DC 

220  N'  = 1.P'RS 

00024000 

M=  P.i-1 

00027000 

MR  r MINO ( IPANH.N R-M) 

00026000 

CP=  RIP!) 

00027000 

R<N)  = CP/AK  <N  . 1 ) 

00020000 

DC 

220  | r 2 . MR 

00027000 

I = M ♦ L 

OOOJCCOO 

220 

R(T)=  PIT)  - AK(N,l)«CP 

00021000 

R ( NIP  ) r R(MR)/AK(NR.l  ) 

00072000 

DO 

320  I : l.NRS 

00077000 

Nr  P'R-  I 

00074000 

Mr  M«  1 

0007*000 

MR  r MINJO  ( TRAP  P .P.R-M  ) 

00076000 

OC 

*20  K r 2. MR 

00077000 

L r M +K 

O0O7POPP 

STORF 

COMPUTER  PTSPLACFMFP.TS  IP.  LOAD  VECTOR  R 

0007*000 

*20 

R ( P' ) r PIN1)-  AK(N,K)*P(L) 

00040000 

400 

RETURN 

0004 i ooo 

ENC 

00042000 

L.l. 


J 


